1 INTRODUCTION 



Insights into the mechanisms of electromediated gene delivery and 
application to the loading of giant vesicles with negatively charged 
macromolecules 1 ^ 
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We present experimental results regarding the electrotransfer of plasmid DNA into phosphatidylcholine giant unilamellar vesicles 
(GUVs). Our observations indicate that a direct entry is the predominant mechanism of electrotransfer. A quantitative analysis 
of the DNA concentration increments inside the GUVs is also performed, and we find that our experimental data are very well 
described by a simple theoretical model in which DNA entry is mostly driven by electrophoresis. Our theoretical framework 
allows for the prediction of the amount of transfered DNA as a function of the electric field parameters, and thus paves the way 
towards a novel method for encapsulating with high efficiency not only DNA, but any negatively charged macromolecule into 
GUVs. 
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1 Introduction 

Lipid membranes are essential constituents of living organ- 
isms, especially because of their impermeability to ions and hy- 
drophilic molecules. This impermeability allows for the crucial 
compartmentalization necessary for life to develop 1 . However, 
the barrier presented by the membrane turns out to be a hurdle 
in many biotechnological applications such as gene delivery or 
encapsulation of charged compounds into giant liposomes, i.e. 
when one needs to force the way of a charged molecule through 
a lipid bilayer. 

As far as the loading of giant unilamellar vesicles (GUVs) is 
concerned, the most common procedure is to directly prepare 
the liposomes in a medium containing the desired compound. 
Fabrication of vesicles trapping salt-containing solutions can 
be achieved e.g. by hydration of a hybrid film of lipids and 
agarose". Also, natural swelling of GUVs encapsulating high 
molecular mass DNA was shown to be possible in presence of 
moderate concentrations of magnesium ions 3 . Unfortunately, 
these techniques still suffer from the drawbacks inherent to the 
gentle hydration method: the resulting vesicles exhibit a wide 
size heterogeneity and their unilamellarity can not be guaran- 



teed. The well known electroformation protocol^ has a higher 



f Electronic Supplementary Information (ESI) available: one movie and asso- 
ciated legend. See DOI: 10.1039/bOOOOOOx/ 

" Institut de Pharmacologic et de Biologic Structurale, CNRS UMR 5089, Uni- 
versite Paul Sabatier, Toulouse, France. 

b Laboratoire de Physique Theorique, CNRS UMR 5152, Universite Paul 
Sabatier, Toulouse, France. 

c Institut Fresnel, CNRS UMR 6133, UniversitCs Aix-Marseille, Marseille, 
France. 

* Tel: +33 5 61 17 58 11; E-mail: rols@ipbs.fr. 



yield of unilamellar vesicles 5 . This protocol was recently 
refined in order to make electroformation possible with solu- 
tions containing various compounds, and in particular signifi- 
cant amounts of ions. Performing the electroformation in a flow 
chamber allows the encapsulation of high ionic strength solu- 
tions or large dextrans in the so formed GUVs 6 . However, this 
technique requires the use of a sophisticated setup and does not 
seem easy to implement. Recently proposed modifications of 
the electroformation protocol enabled the formation of GUVs 
under physiological conditions^. Nevertheless, these methods 
have not been tested with large compounds of several tens of 
kDa, nor with highly charged molecules. Another solution was 
proposed by Stachowiak et ai, who designed a GUV formation 
procedure using a pulsed microfluidic jet that deforms a planar 
bilayer into a vesicle 9 . This technique does produce unilamel- 
lar vesicles of homogeneous radii, and any solution could in 
principle be trapped in the GUVs. Nevertheless, implementing 
the microfluidic setup is not a trivial task. 

For the purpose of gene delivery to living cells, viral based 
methods are available; these are quite efficient but their safety 
has been questioned l(> . Chemical methods relying on the for- 
mation of DNA complexes with positively charged molecules 
are also a topic of active research 1 1 ; such methods, although 
safer, are less efficient than the viral methods. Among physi- 
cal methods, electropermeabilization is one of the most widely 
used. Indeed, the cell membrane can be safely and transiently 
permeabilized in a very elegant manner by applying electric 
pulses 12 . Providing that the pulses are of sufficient durational 
and amplitude 14 , otherwise non permeant molecules of ther- 
apeutic interest can enter the cytoplasm of mammalian cells. 
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This fact has led to two clinical applications: electrochemother- 
and electrogenetherapy^. The former involves small 
molecules such as cisplatin or bleomycin, and the latter larger 
molecules such as plasmid DNA. 

As was described in^, the mechanisms of molecular uptake 
under electropermeabilization, although poorly understood for 
the moment, depend strongly on whether the transfered objects 
have molecular or macro-molecular sizes. Whereas smaller 
ones seem to be able to freely cross electropermeabilized cell 
membranes, macromolecules exhibit an intermediate interac- 
tion with the membrane, the degree of this interaction being 
correlated with the ultimate transfection efficacy. A precise de- 
scription of the electropermeabilized membrane would help to 
design safer and more efficient protocols; this is a major mo- 
tivation for investigating the behaviour of simpler model sys- 
tems. To our knowledge, DNA electrotransfer into liposomes 
was first studied by Chernomordik et alJ^. They found that 
high molecular mass DNA could enter DPPC/cholesterol (7:3, 
mokmol) LUVs, via endocytosis-like vesicles which shielded 
the electrotransferred DNA from the internal medium. The in- 
ternalization mechanism they proposed was based on their ob- 
servation of the characteristic fluorescence of DNA/ethidium 
bromide (EB) complexes after sonication of EB-loaded vesi- 
cles which were pulsed in the presence of DNA. However, they 
could not directly observe the liposomes because of their small 
size. The conclusions of 18 were subsequently questionned by 
Lurquin and Athanasiou 19 , who observed that giant (« 10 /urn) 
EB-loaded DPPC liposomes pulsed with DNA did actually 
show the bright fluorescence of DNA/EB complexes before (or 
even without) sonication. These results support a mechanism 
involving electropores and a direct entrance of DNA into the li- 
posomes. No endocytosis-like vesicles occur during this mech- 
anism, and this direct entrance allows the DNA to be imme- 
diately in contact with the medium inside the liposomes. Al- 
though electropores could not be unequivocally observed, these 
results are at variance with the findings 

In this paper we describe high molecular mass DNA elec- 
trotransfer experiments performed on egg phosphatidylcholine 
(EggPC) GUVs. We present qualitative observations regarding 
the pathway of electromediated DNA entry into the liposomes, 
and quantitative results concerning the evolution of the inner 
DNA concentration as a function of the number, amplitude and 
duration of the applied electric pulses. By comparing our ex- 
perimental data to a simple theoretical model, we are able to 
address the two following key questions: (i) what is the mecha- 
nism of DNA electrotransfer across pure lipid membranes? and 
( ii) is it possible to control/predict the amount of electrotrans- 
fered DNA? As well as being of interest for its comparison with 
the uptake of DNA by living cells, this study presents a method 
for loading vesicles with negatively charged macromolecules 
which may have a number of academic and practical applica- 
tions. 



2 Materials and methods 

2.1 GUVs 

Egg yolk L-a-phosphatidylcholine (EggPC) and L-oc- 
phosphatidylethanolamine-N-(lissamine rhodamine B sul- 
fonyl) (Rhodamine PE) were purchased from Avanti Polar 
Lipids (Alabaster, AL). Lipids were diluted in chloroform, 
at a mass concentration of 0.5 mg/mL, and stored at -20 ° C. 
Rhodamine PE dye was added at a concentration of 1 mol %. 

The vesicles were prepared at room temperature using the 
electroformation protocol 4 , in an aqueous solution of 240 mM 
sucrose (internal solution) and subsequently diluted in an aque- 
ous solution of 260 mM glucose and 1 mM sodium chloride 
(external solution). The pH of the external solution was ad- 
justed to 7 using a 1 mM phosphate buffer (KH2PO4/K2HPO4, 
Merck, Darmstadt, Germany). The pH of the internal solu- 
tion was measured to be 6.6. This slight initial pH asymme- 
try did not have any effect on our experiments: non-pulsed 
vesicles were stable for several tens of minutes, and no in- 
crease in DNA concentration inside a non-pulsed GUV could 
be detected. Furthermore, this initial pH difference will van- 
ish after a few permeabilizing pulses because of the mixing of 
inner and outer solutions. The conductivities of the internal 
and external solutions were measured with conductivitymeter 
HI 8820 (Hanna Instruments, Lingolsheim, France), and had 
values of approximately 1.5 x 10~ 3 S/m and 4.5 x 10~ 2 S/m, 
respectively. The osmolarities, measured with Osmomat 030 
osmometer (Gonotec, Berlin, Germany), were approximately 
280 mOsm/kg and 300 mOsm/kg, respectively. These solu- 
tions are typically used when working with GUVs for several 
reasons. The similar osmolarities ensure that the GUVs have 
very little initial membrane tension, and the low conductivi- 
ties prevent Joule heating in the pulsation chamber. The sugar 
asymmetry yields a density difference that allows sedimenta- 
tion of the vesicles to the bottom of the chamber and facilitates 
their localization. 

Electroformation was performed as follows. A small volume 
(15 /jL) of the lipid solution in chloroform was deposited on the 
conducting sides of glass slides coated with indium tin oxide. 
The glasses were then kept for two hours under vacuum in a 
desiccator to remove all traces of the organic solvent. After- 
wards, the plates, spaced by a 1 mm thick silicon frame (Elec- 
tron Microscopy Sciences, Hatfield, PA) were assembled into 
a chamber. The chamber was filled with the sucrose solution 
(internal medium). The slides were connected to an AC field 
function generator (AC Exact, model 128; Hillsboro, OR) and 
sinusoidal voltage of 25 mV peak to peak at 10 Hz was applied. 
The voltage was increased by 100 mV steps every 5 minutes, up 
to a value of 1225 mV and maintained under these conditions 
overnight. Finally, square-wave AC field of the same amplitude 
was applied at 5 Hz for one hour in order to detach the GUVs 
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2.5 Image and data processing 



from the slides. 

2.2 DNA 

A 4.7 kbp plasmid (MW - 3.1 6 Da) pEGFP-Cl (Clontech, 
Mountain View, CA) carrying the green fluorescent protein 
gene controlled by the cytomegalovirus promoter was stained 
stoichiometrically with the DNA intercalating dye TOTO-1 
(Molecular Probes, Eugene, OR). Staining was carried out at 
a DNA concentration of 1 /Jg//jL for 60 minutes on ice, at a 
basepair to dye ratio of 20. Plasmids were prepared from Es- 
cherichia coli transfected cells by using Maxiprep DNA purifi- 
cation system (Qiagen, Chatsworth, CA). The plasmid was then 
diluted in the pulsation buffer at a mass concentration of 1 /jg 
per 100 /jL, a value typically used in cell electrotransfection 
experiments. This plasmid was chosen because of its large size 
similar to those of the DNAs used inl 18 l 19 [ and also because it 
is classically used in gene electrotransfer studies where it was 
observed to form discrete interaction sites at the cell surface 
during electropermeabilization experiments^! 

2.3 Electropulsation 

Pulsation chambers were similar to those used in 20 . Two paral- 
lel copper strips (3M, Cergy-Pontoise, France) were stuck on a 
glass slide 0.5 cm apart. A glass coverslip was then stuck onto 
the glass slide with heated parafilm. The cavity between the 
slide and the coverslip was filled with 30 /jL of the buffered glu- 
cose solution and 2 /jL of the GUV solution. Electropulsation 
was performed directly on the microscope stage, by applying a 
series of 10, 20 or 40 pulses of 0.5, 1 or 5 ms duration at 0.33 
Hz repetition frequency with a p-Tech pulse generator ((3-Tech, 
L'Union, France). Field amplitudes Eq were adjusted according 
to the initial vesicle radius Rq in order to impose initial induced 
transmembrane voltages A\\Iq = (3/2)/?o^o 21 ranging between 
0.5 and 2 V, and as a consequence varied between 20 and 80 
kV/m from one experiment to another. 

2.4 Confocal microscopy 

All experiments were performed under an inverted confocal mi- 
croscope (Zeiss LSM510; Carl Zeiss, Jena, Germany) equipped 
with a 63 x Zeiss objective for fluorescence imaging. For the 
red channel (membrane labeled with Rhodamine), excitation 
at 543 nm was provided by a HeNe laser, and emission fil- 
ter was a 560 nm long pass. For the green channel (DNA la- 
beled with TOTO-1), excitation at 488 nm was provided by a 
Ar laser, and emission filter was a 500-530 nm band pass. Im- 
ages were acquired sequentially. Acquisition time was 2 s (1 s 
for each channel), thus one image on each channel could be ac- 
quired between every two consecutive pulses, separated by 3 s. 
The characteristic time for homogenizing the concentration of 
a molecule of diffusion coefficient D within an object of size 



/ is on the order of l 2 /D. For plasmid DNA of diffusion co- 
efficient D = 10~ 12 m 2 s -1 l22land regions of inhomogeneous 
DNA concentration with sizes in the micrometer range, this 
yields characteristic times on the order of a few seconds. Thus 
in our experiments, DNA concentrations had enough time to 
homogenize between consecutive pulses. The use of confocal 
microscopy for these experiments is crucial, as it ensures that 
the measured fluorescence intensities are simply proportional 
to the local concentration of the fluorescent molecules 23 . 

2.5 Image and data processing 

We were interested in two different quantities, both measur- 
able via confocal fluorescence microscopy: the vesicle radius 
R and the ratio of the DNA concentrations inside and out- 
side the vesicle, c/cq. As we work in a dilute water solutions 
with similar physicochemical properties, the concentration of 
a fluorescently-labeled molecule is simply proportional to the 
mean fluorescence intensity due to this molecule, thus the quan- 
tity c/cq is equal to I /Iq, the ratio of the mean fluorescence 
intensities inside and outside the vesicle, which is easily mea- 
surable in our experiments. 

Our images were composed of two channels: a red one for 
Rhodamine PE (lipid membrane) and a green one for TOTO- 1 
(DNA). GUVs size measurements were performed on the red 
channel in a semi-automatic manner via custom-written pro- 
grams, and relative DNA amount quantification was performed 
on the green channel automatically, by computing the mean 
fluorescence intensity inside the vesicle within a disk whose 
radius was chosen at a value of 80 % of the measured radius 
of the vesicle (we checked that no significant differences were 
seen for disk size values ranging from 60 to 90 % of the GUV 
size). This procedure was repeated for each image, so we were 
able to plot the vesicle size and the DNA quantity trapped in 
the vesicle as a function of time. We also computed the mean 
fluorescence intensity outside the vesicle. It was checked that 
this value did not change significantly after each pulse, and thus 
we used the value computed from the first image. 

All images had non-zero mean fluorescence intensity values, 
even inside a non-pulsed vesicle where no DNA was present, or 
in pictures of pulsation chambers filled with pure water. This 
was due to the bias level of the confocal microscope sensors. 
In order to properly quantify DNA concentrations, we esti- 
mated this bias fluorescence level inside a non-permeabilized 
GUV containing no DNA, and substracted this quantity from 
the measured internal and external DNA fluorescence intensi- 
ties to obtain the fluorescence intensity ratio I/Iq equal to the 
concentration ratio c/cq. 

Data for c/cq exhibited rather strong fluctuations and were 
smoothed using a n = 3 moving average technique. Error 
bars in Fig. [2] were obtained by computing the standard devi- 
ation of the fluctuations around the constant value of c/cq af- 
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ter the application of pulses, for each experiment. Image and 
data processing tasks were performed with Matlab (The Math- 
works, Natick, MA) and ImageJ (National Institute of Health, 
Bethesda, MD). 

2.6 Numerical computation of the electric field near a per- 
meabilized vesicle 

As described in the theoretical section of this paper, we were 
interested in flow of DNA into a permeabilized vesicle. In 
the limit where the electrophoretic force on a DNA molecule 
dominates the osmotic force (due to concentration gradients), 
which we argue is the case in our experiments, the entry of 
DNA into the cell can be computed from the flux of the elec- 
tric field across the cathode-facing permeabilized region of the 
GUV. The permeabilized region is defined as a spherical cap of 
angle facing the cathode, and we also assume that an identical 
region, similarly permeabilized to ionic currents, exists on the 
anode-facing side. 

We computed numerically this flux (]) as a function of 8 us- 
ing the finite element calculus software Comsol Multiphysics 
(Comsol, Burlington, MA). 

The GUV was modeled as a spherical shell of radius R and 
internal electrical conductivity a,, embedded in an aqueous so- 
lution of electrical conductivity a e . Laplace's equation for the 
electric potential 



V • aV\|/ = 



(1) 



was solved numerically in a cylindrical box of radius L and 
heigth 2L, taking advantage of the axial symmetry of this prob- 
lem. An electric field of amplitude Eq = Vq /L was applied in 
the z-direction parallel to the cylinder axis, by imposing the 
Dirichlet boundary conditions l|/ = Vq and \|/ = — Vq on the top 
and the bottom faces of the cylinder, hence their names of an- 
ode and cathode, respectively. The membrane was modeled 
using thin layer boundary condition with a thickness d = 4 nm 
and a membrane conductivity a m = 5 x 10~ 7 S/m 25 . For repre- 
senting the permeabilized regions of the vesicle facing the cath- 
ode and anode, we replaced the membrane conductivity o,„ by 
the external conductivity o e on a spherical cap of angle 9. Nu- 
merical integration of E = — V\|/ on that surface thus gave us 
the value of the flux (j) as a function of 9. 

Different values of the model parameters R, Eq, o, and 
G e corresponding to realistic experimental conditions were 
explored, and are given in the theoretical section and in 
Figs.|]and|4] Typical values were R = 10 5 m, Eq — 50 kV/m, 
Oi = IQ- 2 S/m and o c =4x 10~ 2 S/m. 

2.7 Determination of the relative concentration incre- 
ments 8x 

The amount of DNA transfered into a GUV per pulse can be 
considered in terms of a relative concentration increment 8x, 



i.e. the difference between the value of the DNA concentration 
inside a liposome after a pulse and its value before the pulse, 
normalized by the outer DNA concentration. This quantity was 
measured experimentally by substracting the value of the mean 
fluorescence intensity ratio I/Iq before the pulse to the value of 
I/Iq after the pulse. 

Within the framework of our model, Sx can also be pre- 
dicted using Eq. 



10 



where the electric pulse parameters Eq 
and x are chosen by the experimenter, the vesicle size R is 
easily measured, and the DNA electrophoretic mobility is /j = 
-3.75 x 10~ 8 m 2 V _1 s -1 E§l The value of the flux factor/' (0)9 
is estimated to be on the order of 0.2 (see section|4]and Fig. [4]). 



3 Experimental results 

We performed DNA electrotransfer experiments with 21 
EggPC liposomes. Either 10, 20 or 40 pulses were applied at 
0.33 Hz. The durations of the applied pulses were 0.5, 1 or 5 
ms, and correspond to those commonly used for electromedi- 
ated gene transferal The longest pulse duration x was observed 
to maximize the amount of transfered DNA, so most of the ex- 
periments were done at this value of 5 ms. Depending on the 
experiment, the pulse amplitude varied from 20 to 80 kV/m. 
The applied field magnitude Eq was tuned in order to obtain 
initial induced transmembrane voltages (at the poles facing the 
cathode and the anode) ranging between 0.5 and 2 V. This was 
estimated from the formula for the voltage drop at a point on 
the membrane whose radial vector (of magnitude R) makes an 
angle 9 with the direction of the applied fieldEB 



A\\f - 



-E Rcos{d). 



(2) 



This formula is valid when the conductivity of the membrane is 
much smaller than those of the internal and external solutions 
and when the thickness of the membrane d is small compared 
to the vesicle size R, which is true here. In some cases we 
were able to keep several vesicles in the optical field, thus we 
could gather data from several liposomes during the same ex- 
periment. Typically, we tuned the field magnitude according to 
the size of the largest GUV we could see. IrJ22l it was estab- 
lished for EggPC and DOPC liposomes that there was a critical 
value of the absolute value of the transmembrane voltage A\|/ c 
below which vesicles were not permeabilized and showed no 
change under the application of the electric field. Sufficiently 
large vesicles were seen to eject lipid material in the form of 
small vesicles and tubules, however this lipid ejection caused 
the vesicle radius to diminish and at constant field, when the 
radius R fell below a threshold value the vesicle was no lon ger 
visibly affected by the field. In agreement with the results or 2 ^ 
it was found that the smallest liposomes were not permeabilized 
and no DNA fluorescence could be detected inside them. For 
vesicles of intermediate sizes, the permeabilized area is smaller 
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than that of the largest GUV on which the field of view was cen- 
tered. Consequently these vesicles of intermediate sizes should 
be expected to exhibit relatively low amounts of DNA uptake. 
However, we kept such vesicles in our analysis, and coherent 
results for the increments of the internal DNA concentration 
were also obtained with them (section[4]). 

As can be seen in Fig. [T] it is possible to optically detect 
DNA entering the vesicle - a movie is also provided as ESIf on 
the journal website. The images show vesicles exhibiting lipid 
loss via formation of tubular structures. These structures were 
created on the anode facing side of the GUV and appeared to 
remain attached to the vesicle and stable over a few minutes, as 
previously reported in 20 for DOPC liposomes. 

Darker zones denoting a local depletion of DNA from the 
outer region of the vesicle facing the cathode were observed in 
a small number of experiments, as shown in Figs.[T]3 and C. 
This phenomenon was quite rare, and could only be detected 
following the first few electric pulses. It may be due to hydro- 
dynamic flows resulting from the leakage of the inner sucrose 
solution or from osmotic effects. However, this local deple- 
tion did not have any appreciable influence on the subsequent 
behavior of the vesicle around which it was seen. 

We confirm the results of Lurquin and Athanasiou 19 , namely 
the free entrance of DNA through pores formed by the elec- 
tric field. Indeed if an uptake via endocytosis-like vesicles was 
the predominant mechanism, we would detect the presence of 
such objects in most of our experiments. Even though presum- 
ably quite small (some tens or hundreds of nm), some of these 
endocytosis-like vesicles should be created in the focal plane 
where the transmembrane voltage induced by the electric field 
is maximal, and should thus appear in our images. This was not 
the case. Therefore our pictures clearly speak in favour of an 
entrance via pores. We mention that we sometimes noticed en- 
docytosis events, but they were very rare, and these phenomena 
always took place after application of the electric pulses. A de- 
tailed explanation of this rare uptake via vesicles is beyond the 
scope of the present work. However, we stress that an explana- 
tion based on the mechanism described in^ where long DNA 
fragments of a few kbp absorb on a patch of the GUV mem- 
brane before it bends and pinches off toward the interior of the 
GUV would not be satisfactory because in that study the au- 
thors reported no DNA internalization within liposomes made 
just of zwitterionic lipids, as our EggPC GUVs. 

Typical examples of fluorescence experimental data, ob- 
tained as described in the Materials and Methods section, are 
shown in Fig. [5] (diamond marks). These data for I/Iq are the 
ratio of the mean fluorescence inside and outside the GUV, and 
are equivalent to the relative DNA concentration c/cq, where c 
is the DNA concentration inside the liposome and cq the DNA 
concentration outside the liposome. A ratio I/Iq equal to 1 thus 
means that the DNA concentrations inside and outside the vesi- 
cle are the same. The field amplitude Eq, the pulse duration 



T, the number of pulses N, the initial vesicle radius Rq and the 
initial transmembrane voltage A\\tQ induced by the electric field 
are given in each panel. For each experiment, we observe an in- 
crease of the amount of DNA inside the vesicle during the train 
of pulses, along with a stabilization after the end of the pulse 
sequence. This clearly shows that DNA can enter GUVs and 
remain trapped in a stable manner, without significant leakage. 
We draw the reader's attention to the different scales used on 
each panel: the final DNA concentrations inside the liposomes 
presented in the two top panels were much lower than those 
inside the liposomes subjected to electric pulses of the longest 
duration, which are presented in the bottom panels. 

A number of control experiments were performed. First, 
we checked that fluorescence intensities inside the liposomes 
did not change if we did not apply pulses. This verification 
was performed on long timescales (several minutes). As a sec- 
ond control experiment, we applied a series of ten pulses and 
checked that the fluorescence level, after having increased, re- 
mained stable until the application of another train consisting of 
ten pulses of slightly higher amplitude, which induced a subse- 
quent augmentation of the relative fluorescence intensity. This 
is important as it proves that the DNA uptake is induced directly 
by the electric field, ruling out the possibility that the vesicle re- 
mains permeable to DNA for some time after the field is cut and 
emphasizing the importance of the electrophoretic force on the 
DNA as well as the permeabilizing effect of the electric field on 
the membrane. We also carried out control experiments with 
TOTO-1 but with no DNA, in order to check that it was not 
unbounded dye we observed. The quantum yield of TOTO-1 is 
dramatically increased when bound to DNA. As a consequence, 
images acquired with same microscope settings and in the ab- 
sence of DNA should be much less fluorescent than those with 
DNA. This was indeed the case (in fact the signal from TOTO- 
1 alone was too weak to be detected), thus proving that it was 
indeed fluorescently labeled DNA we were monitoring. 

Interestingly, we found that sometimes the final DNA con- 
centration inside the GUV could exceed the outer DNA con- 
centration (c/cq = I/Iq > 1). That can be seen for example on 
the bottom right panel of Fig. [2] This is quite surprising, and 
means that if the DNA enters the vesicle from the cathode fac- 
ing pole driven by the electric field, it is not able to freely leak 
out the vesicle from the other pole. Indeed, if the permeabi- 
lization structures were the same on each hemisphere of the li- 
posome, the maximal DNA concentration ratio possible should 
be c/cq = 1. This is not the case in our experiments. This 
observation, along with the asymmetric formation of tubules, 
highlights the different membrane reorganization and symme- 
try breaking between the two hemispheres. The possibility of 
exceeding the external concentration of DNA within the vesi- 
cle also implies that the electrophoretic force due to the ap- 
plied field overcomes the osmotic force due to the difference in 
the concentrations between the inside and outside of the vesi- 
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Fig. 1 A-F: Typical raw images of one experiment. Rhodamine PE (lipid membrane) and TOTO-1 (DNA) are shown in red and green 
respectively. 20 pulses of 39 kV/m amplitude and 5 ms duration were applied at 0.33 Hz (one pulse every three seconds, and a total pulsation 
time of 60 s). Images were taken at 0, 6, 12, 18, 54 and 222 s; time origin was defined as the onset of the acquisition, and the first pulse was 
applied after about 2.5 s. Vesicles are thus shown after having received 0, 2, 4, 6, 18 and 20 electric pulses, respectively. Field polarity is 
indicated on pictures acquired during electropulsation (+ and - electrodes). The last picture (panel F) was acquired more than two minutes 
after the pulse train, and shows the stability of the vesicles loaded with DNA. Tubular structures associated with vesicle size decrease and 
similar to those described in 20 are present. A scalebar of 10 /jm length can be seen on panel F. G and H: Fluorescence intensity profiles along 
white rectangular regions shown in panels A and F, respectively. Thick green line represents the mean value of TOTO-1 fluorescence intensity 
inside, on the left, and on the right of the GUV. Thick black line represents the bias fluorescence level, a non zero fluorescence value 
corresponding in reality to a zero DNA concentration. This value was thus substracted for subsequent data analysis. 



cle, this will be discussed in the theory section which follows. 
We were also able to monitor the size decrease of the GUVs 
during the experiments. Applying the method previously de- 
scribed in^ we obtained the fraction of permeabilized area 
lost per pulse X and the critical transmembrane potential differ- 
ence A\\t c required to observe the size decrease of the GUVs. 
For the pulses of 5 ms duration (the same duration x as used 
irJ^SI), we find by averaging on all our experiments: X ps 0.32 
and A\\t c pa 0.50 V. The value found for X agrees well with the 
previously reported value of 0.31, which means that the GUVs 
still lose ~ 30 % of their permeabilized area (area where A\|/ c is 
exceeded) per pulse 20 . However, the A\|/ c value of 0.50 V found 
here is much smaller than the value 0.89 V found for vesicles 
pulsed in the absence of DNA. This suggests that the presence 
of DNA in the external medium, and/or inside the vesicle af- 
ter a few permeabilizing pulses, tends to destabilize the vesicle 
membrane and makes it easier to electropermeabilize, but does 
not affect the fraction of permeabilized area expelled from the 
liposome. 

4 Theoretical analysis of the relative DNA con- 
centration increments 

The theoretical interpretation of our experimental results is 
rather difficult as we do not have explicit access to the way 



in which the vesicles are permeabilized. 

In this section we will present a simple model describing the 
relative DNA concentration increment 8x = 6c/cq inside a vesi- 
cle of radius R, caused by an electric pulse of amplitude Eq and 
duration x. The DNA concentration increment inside the vesi- 
cle, per pulse, is denoted by 8c, and the constant outside con- 
centration by co. We will assume that the major contribution 
to 8x is due to the pulling of the DNA molecules by the elec- 
tric field, and we will see that our 360 experimental data points 
agree with this hypothesis. 

The DNA concentration c is governed by an electrodiffusion 
equation 



where j is the thermodynamic current j = — DVc + /jcE. The 
first term is a diffusion term that depends on the molecule dif- 
fusion coefficient D. The second term is an electrophoretic term 
depending on the local field E and the electrophoretic mobility 
of the molecule /j (the steady state velocity v of the molecule 
in a field E is given by v = fjE). The DNA concentration incre- 
ment 8c(f ) inside an object 1/ of volume V during a time t is 
given by 
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E Q = 67 kV/m x = 0.5 ms N = 20 pulses E Q = 28 kV/m x = 1 ms N = 40 pulses 

R Q = 19 urn At|f Q = 1 .93 V R Q = 22 urn A\|/ Q = 0.91 V 




time [s] time [s] 

Fig. 2 Plots of the relative mean fluorescence intensity inside the vesicle I/Iq (or equivalently the relative concentration of DNA c/cq) as a 
function of time, for 4 typical experiments. Diamond marks represent raw data, and thin full lines these same data smoothed by a n — 3 
moving average procedure. This smoothed values were used to compute the individual concentration increments &c. Error bars were obtained 
by computing the standard deviation for the I/Iq values after the pulse sequence. Final relative concentration, obtained by averaging the I/Iq 
values after the pulse sequence, is shown with horizontal dashed lines. Thick vertical lines indicate the end of the pulse sequence. Note the 
different horizontal and vertical scales for each graph. Electric field parameters, initial vesicle radius Rq and initial induced transmembrane 
voltage A\|/q are given on each panel. One pulse was applied between every two consecutive data points (0.33 Hz repetition frequency). 
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which can be written using Gauss' theorem: 




We easily estimate the relative importance of electrophoretic 
forces and DNA diffusion by comparing the diffusive or os- 
motic part of the current per concentration r os = j os /c = 
— DVln(c) and the electrophoretic part r ep =j ep /c = /jE. We 
use the following values of DNA in solution: /j = —3.75 x 
10" 8 m 2 V" 1 s" 126 andD= l(T 12 mV 122 . We use the estima- 
tions |r M | ~ D\ln(c/co)\/R (where c is the average DNA con- 
centration in the vesicle and R the vesicle radius) and \r ep \ ~ 
\/j\Eo, this gives 

|r J D\\n(c/c )\ 4 fc\ 

\r ep \ \n\E R \cqJ 

for the lower range of permeabilizing applied voltages (EoR ~ 
1 V, cf subsection |2.6| l. We thus see that the effects of diffusion 
are negligible except for perhaps the first pulse when c ~ 0. 

As far as the ionic currents are concerned we assume that the 
vesicle is permeablized in a symmetric way at both the cathode 
and anode facing sides. The permeabilized region subtends an 
angle in the direction of the applied field. As mentioned in 
the Material and methods section we model the permeabilized 
region by replacing the membrane conductivity by that of the 
external solution; this is clearly a simplification. As far as the 
transport of DNA is concerned the fact that macropores are only 
ever seen at the cathode facing pole suggests that although the 
conductive properties of the permeabilized regions at the an- 
ode and cathode facing poles are the same (i.e. the transport 
properties of these regions to small ions is the same), DNA can 
pass more freely at the cathode facing side than the anode fac- 
ing side. Thus in our analysis we posit that DNA can be freely 
electrophoretically conducted through the cathode facing per- 
meabilized region. In our model this is equivalent to saying 
that electrophoretic mobility of DNA in the membrane region 
is zero except in the permeabilized region facing the cathode, 
where it takes its free solution value. The concentration incre- 
ment 8c (f ) is simply related to the integral over time of the flux 
of the thermodynamic current j across the surface S. Let us 
now consider a spherical object with fixed volume V and let us 
compute the concentration increment 8c from the onset of one 
pulse to the onset of the subsequent pulse. Using a constant 
volume is a reasonable approximation because in our experi- 
ments the average volume decrease caused by a single pulse is 
still small compared to the initial volume of the vesicle. 

We will base the analysis of our experimental results on the 
following assumptions (a) the dominant transport process is 
electrophoretic and we neglect DNA diffusion (b) we assume 
that over the period of the pulse a conducting hole is created 
at both poles of membrane each having radius a. Concretely 



this means that the average radius of the conducting hole in the 
membrane averaged over the time of the pulse application is 
a. Beyond the time averaging, this is a crude simplification, as 
the surface permeable to ions may be larger than the one per- 
meable to DNA, this is because we cannot be sure that there is 
really a single macropore in the lipid bilayer, instead of several 
smaller pores, or other types of defects. Nevertheless, we ex- 
pect the electroconvected flux of DNA into the vesicle to be of 
the same order. Furthermore we will assume that DNA is not 
able to exit the vesicle at the conducting region facing the cath- 
ode. Adopting the hypotheses above, upon a pulse of duration 
x the increase in concentration is given by 

5c= \j^f E . dS . ( 7) 

i.e. concentration increment 8c is just proportional to the flux 
(j) of the local electric field E across the permeabilized surface 
S p , the hole of radius a. In terms of the relative concentration 
increments 8x = 8c/co, which are exactly the differences be- 
tween every two consecutive points of our experimental data 
for I/Iq, this can be written 

S*=»4>. (8) 
4%R 3 Y 

In order to be able to evaluate 5x, we should be able to evaluate 
(j). Recall we assume that the permeabilized area is a spheri- 
cal cap of angle facing the cathode, that both DNA and ions 
can cross. Despite this simplification of the transport prob- 
lem no analytical solution to this problem exists. However, on 
purely dimensional grounds (assuming the limit of zero mem- 
brane thickness), we can write 

4 = E Q R 2 f(a/R) , (9) 

where a is the pore radius such that a = Rsm(Q), and / an 
unknown function of the ratio a/R, or in other words a function 
of the permeabilized angle 0. Note that in the limit where the 
pore length is much greater than the pore radius, which is the 
case for nanopores, Eq. [9] will not be valid 28 and will depend 
explicitly on the pore length. To check the validity of Eq. [9] 
we computed numerically (j) for different field intensities Eq, 
different vesicles radii R, and different angles 0, as described in 
the Materials and methods section. These results are presented 
in Fig. [3] where one can see that all the flux values superpose. 
These plots are in fact plots of the function f(a/R) = /(0), and 
show that Eq. [9] is valid. For a small angle 0, we can make 
the linear approximation that /(0) ps /'(O)0, where /'(0) is the 
slope of the straight line giving /(0). We finally obtain for the 
relative concentration increments 

&=»/(0)9. (10) 



8 



4 THEORETICAL ANALYSIS OF THE RELATIVE DNA CONCENTRATION INCREMENTS 
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Fig. 3 Flux (]) of the electric field across a permeabilized spherical 
cap of angle 6. Flux values were computed numerically for different 
vesicle radii R ranging from 5 to 50 /jm, and various electric field 
amplitudes Eq ranging from 10 to 500 kV/m. Internal and external 
conductivities <T, and o e were both equal to 4 x 10~ 2 S/m. Slight 
deviations occur for curves corresponding to the different radii, and 
are due to the different meshes used. All curves for the different field 
amplitudes exactly superpose. 

We checked the approximation /(8) « /'(O)0 by fitting 
straight lines to the numerically computed /(8). As shows 
Fig. [4] the fits are in very good agreement with numerical data. 
The different plots correspond to different internal to external 
conductivitiy ratios <7,/<7 e , as indicated in the legend. One can 
see that the slopes /'(0) depend on these conductivity ratios, 
and increase for increasing ratios. This makes sense, as one 
would not expect any flux of the electric field for a perfectly 
non conductive object. This is also interesting because it gives 
another reason why the relative DNA concentration increments 
seems to increase during one experiment. This is because the 
solution inside the GUV being initially less conductive than 
the external medium, the electropermeabilization leads to the 
mixing of the internal and external solutions, and thus to an 
increase of a, towards o e . This increase of the conductivity 
ratio causes in turn an increase of the flux (j), and hence an 
increase of Sx. One can also see on Fig. |4] that for conduc- 
tivity ratio values between 1 and 10 corresponding to our ex- 
periments, the flux factor /'(0)8 should be, for moderate an- 
gles, on the order of 0.1-0.5. Its value should of course be 
varying from one data point to another, mostly because of the 
difference between the permeabilized areas, but our analysis 
nevertheless states that if our assumptions are valid, the values 
of /'(0)8 = 4nRdx / (3\/j\Eqz) obtained with our experimental 
data should be distributed around 0.1-0.5. We have gathered 
in Fig. [5] all our experimental values of /'(0)8, using the value 
of /j given above. This histogram contains more than 300 in- 
dividual data points, corresponding to various vesicle sizes R, 
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Fig. 4 Flux <j) of the electric field across a permeabilized spherical 
cap of angle 0. Flux values (dashed lines) were computed 
numerically for different internal conductivities o, of values a e , 
a e /2, a e /3 ... O e /10, G e /20, fj e /40 ... O e /100 and o c /1000 (from 
top curve to bottom curve). Vesicle radius R was 10 /jm, field 
amplitude Eq was 60 kV/m, and external conductivity a e was 
4 x 10 -2 S/m. Thick full lines are linear fits to numerical data. These 
fits show that the linear approximation on /(0) made to obtain Eq.[To| 
is valid, and that for reasonable angles of the permeabilized area 6 
and for internal and external conductivities of the same order, /'(0)6 
is of the order of 0.1-0.5. 



field amplitudes Eq, and pulse durations T. The narrow shape 
of the distribution around the expected values thus supports our 
hypothesis that DNA entry in the GUVs mostly happens via 
electrophoretic effects. 

The value of the flux factor /'(0)6 averaged on all experi- 
ments is (/'(0)8) « 0.26. If we take the value of 2.34 rad -1 for 
/'(0), corresponding to a conductivity ratio equal to 1 (this has 
to be true after a few permeabilizing pulses, once the internal 
and external solutions have mixed), we obtain for the average 
value of the permeabilization angle 9: (8) « 6.3°. This angle 
value seems reasonable, but should not be interpreted as the 
most frequent permeabilization angle in our experiments. In- 
deed, this quantity is extracted from an averaging process on 
many experiments performed with various electrical parame- 
ters and on GUVs of different sizes. Given the fact that the 
vesicles shrank during an experiment, one would expect that 
6 decreases too. Let us recall that in similar systems, single 
macropores could not be detected during the pulse applica- 
tion! 2 ^. The result of this analysis is that the experimentally 
observed DNA transfer into vesicles can be interpreted as, on 
average, that which would be given by electrophoresis through 
an aqueous pore subtending an average angle of about 6° in 
the direction of the applied field. We see that this conclusion 
has the merit of being consistent with the assumptions of our 
model, the effective region across which DNA is transfered is 
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Fig. 5 Histogram of the values of the flux factor /'(O)0 obtained 
from experiments. Due to fluctuations in measurements and to the 
difficulty to detect very low concentration increments at the 
beginning of the pulse sequence, 20 out of 360 experimental 8x 
values were negative and thus lead to negative /'(0)9. Those values 
were removed from this histogram. 



macroscopic but does not represent a large fraction of the vesi- 
cle surface. 



5 Discussion 

Before answering the questions (i) and (ii) introduced in sec- 
tion[T] about the mechanism of DNA electrotransfer into GUVs 
and the ability to control the amounts transfered, respectively, 
we shall briefly discuss our results regarding the decrease in 
size of the GUVs during electropulsation. As observed in!2°J 
the application of electropulses leads to a reduction in the size 
of the vesicles. This reduction in size eventually stops when 
the vesicles have a size such that the permeabilization thresh- 
old is no longer exceeded. Our estimation of the critical mem- 
brane permeabilization potential A\\i c necessary to induce per- 
meabilization turns out to be smaller than that estimated for the 
same vesicle system (EggPC GUVs) in the absence of DNA. 
This suggests the possibility that negatively charged macro- 
molecules may help to induce the permeabilization process via 
some interaction with the membrane, which was also concluded 
from chemical relaxation spectrometry experiments on small 
unilamellar vesicles subjected to an electric pulse in presence 
of DNA 30 . It is possible that DNA may flow into the defects 
that eventually lead to macropores and prevent their closing. 
The presence of DNA in blocked defects could be related to the 
endocytosis based electrotransfer mechanism seen in^. 

(i) Our qualitative observations and the quantitative agree- 
ment of our theory with our experimental results support a 
direct uptake mechanism, as described which is mainly 



driven by electrophoretic effects. However there are some dif- 
ferences between our system and protocol and those used by 
Chernomordik et al. ^ which should be highlighted. For exam- 
ple, in 18 just one pulse was applied whereas we used a train of 
at least ten pulses. As can be seen in the top left panel of Fig. [2] 
the amount of DNA transferred after a single pulse could some- 
times be very small. We could infer that in fact no DNA pene- 
trated the vesicle during the pulse and that the endocytosis-like 
mechanism, as we observed in some rare cases, occured after 
the electrical treatment. The intensities of the field applied 
induced transmembrane voltages of the order of at most 450 
mV whereas ours, and those used by Lurquin and Athanasiou, 
caused initial potential drops of more than 1 V. This could be 
an explanation for why so little DNA has been taken up by the 
liposomes in the experiments of^. Another origin for the dif- 
ference of the reported behaviours could lie in the nature of 
the systems themselves. The liposomes used in-Sl were much 
smaller (LUVs of « 500 nm diameter) and contained 70 % 
DPPC and 30 % cholesterol, whereas the DPPC GUVs o0\ 
and ours made of EggPC, did not contain cholesterol. DNA 
interaction with the membrane of CHO cells during electro- 
transfection experiments appears to be more complex than that 
inferred here. 

InEl 

the initial interaction between DNA and the 
CHO cell membrane leads to the formation of spots on the cell 
surface, again facing the cathode, where DNA aggregates are 
formed. The fact that this spot formation is not seen in GUVs 
has many possible explanations. For DNA and cells it is possi- 
ble that there is a physical or chemical interaction which tends 
to trap DNA near the surface or it could be that the reduced 
mobility of DNA in the cell interior causes a traffic jam like 
phenomenon where the DNA is locally blockedHl. It is worth 
mentioning that the asymmetric transport phenomena discussed 
below also occurs in real The asymmetric transport pat- 

tern found by the authors led them to suggest that pores were 
created on both sides of the membrane, but with a larger size 
(and also lower number) on the cathode-facing hemisphere. 

(ii) The agreement of our experimental results with our 
model for DNA uptake corroborates our quantitative observa- 
tions. It means that we can indeed control and predict the 
quantity of electrotransfered DNA per pulse 8x. Even if we 
can not accurately monitor the size, shape, number and na- 
ture of the membrane defects allowing DNA entrance, we have 
seen that assuming an entry through a macropore of angle 
r~> 6° (or equivalently a flux factor /'(0)9 in Eq. 10 on the 
order of 0.2) provides a satisfactory prediction of the trans- 
fered amounts. Although no macropore could be detected dur- 
ing the pulse application 29 , this simplification turns out to be 
a useful and convenient assumption for estimating 8x. Fur- 
thermore, as was observed in living cells in vitra 33 ' 34 ' and in 
skeletal muscle in viv DNA uptake by the GUVs indeed 
appears to be dominated by the integrated effect of the elec- 
trophoretic force resulting from the electric field application, 
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i.e. the term I^Er/ 1 m Eq. 10 The effect of the electric pulses is 
thus twofold: they both permeabilize the lipid bilayer, and push 
the charged compounds inside the liposome, making the elec- 
trotransfer of highly negatively charged macromolecules very 
efficient. This importance of the electrophoretic effect is rem- 
iniscent of the observation that electric fields generated by salt 
gradients across artificial nanopores enhance the capture rate 
of DNA molecules into the pores^. Another significant exper- 
imental finding is the following. As can be seen on the bottom 
right panel of Fig. [2] we could sometimes reach DNA concen- 
trations inside the vesicle greater than the concentration outside 
(c/co = I/Iq > 1). The inverse dependence of 8x on the vesicle 
radius R alone cannot be responsible for a concentration build 
up within vesicles exceeding that of the bulk. Indeed, even 
if R decreases during an experiment (and thus 8x increases, 
see Eq. [TOj, one should not in principle be able to exceed the 
outer concentration cq. The build up of an excess concentra- 
tion strongly suggests that DNA can not cross the membrane 
on the anode facing pole, and means t hat th e permeant struc- 
tures reported to be created on this sid e 32 ! 36 ! have sizes on the 
order of a few nm at most. It also means that the pore cre- 
ated on the other side closed relatively quickly at the end of 
the pulse, otherwise the inner concentration would always have 
time to equilibrate with the outer one (as the osmotic gradient 
will dominate in the absence of an applied field). Indeed, it 
had been reported that electropores induced by DC pulses in 
GUVs resealed within some tens^l at most a few hundreds of 
ms^ for similar electric pulses, fast enough to prevent signif- 
icant DNA concentration changes due to diffusion. The pres- 
ence on the anode facing pole of the membrane tubules reported 
previously inland which are also observed with EggPC vesi- 
cles (see Fig.[T]i could presumably also prevent the DNA from 
leaving the vesicle from that pole. An interesting corollary of 
these results would be that electromediated loading of vesicles 
with positively charged macromolecules should be much less 
efficient than that observed here. This is because positively 
charged macromolecules would be forced away from the side 
of the vesicles containing the macropores. 

Besides clarifying the mechanism of DNA uptake by vesi- 
cles, our work demonstrates the efficiency of the loading 
method. These results should be relevant to the encapsulation 
of plasmid DNA into giant liposomes, for the purpose of gene 
transfection for example. A simple rule of thumb for choos- 
ing the electric field parameters for a loading protocol would 
be: use a pulse duration of x = 5 ms, tune the field amplitude 
Eq according to the vesicle initial size Rq in order to obtain in- 
duced transmembrane voltages A\|/o = (3/2)RqEq ~ 1 V, and 
then apply a sequence of pulses at 0.33 Hz repetition frequency 
until the required DNA concentration inside the vesicle is at- 
tained. This loading technique is also applicable to other neg- 
atively charged macromolecules, and an equivalent efficiency 
(c/cq > 1) could be attained provided the molecules are big 



enough to be unable to cross the membrane at the anode-facing 
side. The higher the electrophoretic mobility, the higher the 
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electrotransfered amount, as can be understood from Eq 
thus this technique would be particularly efficient with highly 
charged molecules. The negative charge is crucial in order to 
reach the vesicle from the cathode-facing side; indeed, posi- 
tively charged molecules would reach the GUV from the anode- 
facing side, which hosts the nucleation of the tubular structures 
shown in Fig. [T[ and whose crossing is much more difficult, as 
shown by the fact that we oberved surconcentrations of DNA 
(c/c > 1). 

6 Conclusion 

In this paper we have described a quantitative method for load- 
ing GUVs with negatively charged macromolecules, and we 
have shown that the predominant pathway of electromediated 
DNA uptake into liposomes is undoubtedly the electrophoretic 
entrance in a free form via defects created on the cathode-facing 
pole of the vesicles. This spectacular symmetry breaking, first 
observed on GUVs iri^\ is what makes electromediated DNA 
uptake by vesicles so efficient. Indeed, in some of our ex- 
periments, we can reach DNA concentrations inside the vesi- 
cles higher than external concentrations. This would not have 
been possible if large pores also opened on the anode-facing 
hemisphere (as DNA would flow out through these pores), and 
clearly shows that even if the vesicle is permeabilized to small 
compounds at the anode side, sufficiently large moelcules can- 
not cross the membrane in this region. From this we can infer 
that the permeant structures created opposite the positive elec- 
trode have sizes of the order of a few nanometers at most. The 
underlying origin of these striking differences of the membrane 
reorganization at the moleculer level are not yet understood, 
and explaining this asymmetry represents an exciting direction 
for further research. 
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